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摘 要 : 撞击 坑 自 动 识 别 是 行星 探测 任务 中 重要 的 研究 热点 之 一 。 由 于 数据 精度 和 算法 
局 限 性 等 问题 ， 以 往 的 方法 往往 针对 识别 尺度 较 大 的 撞击 坑 。 基 于 目前 国际 上 精度 最 高 、 履 
盖 最 全 的 嫦娥 二 号 数字 高 程 模型 (DEM) 数据 ， 提 出 一 种 基于 等 高 线 自动 识别 小 型 撞击 坑 
的 方法 。 首 先 ， 分 析 小 型 撞击 坑 模 型 的 等 高 线 特 征 ， 确 定 特征 参数 ， 建 立 指标 体系 ; HK, 
提取 实验 区 等 高 线 信息 ， 短 选 出 符合 特征 参数 的 等 高 线 ; 最 后 ， 利 用 霍 夫 变换 拟 合 撞击 坑 边 
界 ， 提取 擅 击 坑 位 置 和 大 小 信息 。 通 过 对 实验 区 域 的 小 型 撞击 坑 识 别 结 果 分 析 表 明 ， 该 方法 
可 针对 不 同 尺度 的 小 型 撞击 坑 设 置 不 同 的 参数 ,分 批 次 识别 提取 ， 综合 识 别 正确 率 达 90% 以 
+, 证 明 该 方法 对 于 自动 识别 提取 小 型 撞击 坑 是 有 效 的 。 该 方法 可 以 弥补 现 有 自动 识别 方法 
往往 针对 尺度 较 大 的 撞击 坑 的 缺陷 ,识别 结果 可 以 补充 现 有 月 球 撞击 坑 数 据 库 中 的 小 尺度 挤 
击 坑 ， 为 今后 建立 全 月 球 撞击 坑 数 据 库 打 下 基础 。 
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行星 表面 通 都 分 布 着 大 量 的 撞击 坑 构 造 , 它们 是 行星 表面 最 显著 的 地 貌 特征 , 也 是 研究 
行星 内 部 物质 的 窗口 中。 由 于 月 球 上 没有 大 气 ， 撞 击 坑 保存 比较 完好 ， 因 此 对 于 研究 月 球形 
貌 、 相 对 地 质 年 龄 、 岩 石 构造 以 及 月 球 的 起 源 有 重要 意义 。 撞 击 坑 识 别 研 究 ， 从 早年 使 用 望 
远 镜 观察 描绘 撞击 坑 的 形 貌 ,到 利用 遥感 图 像 数据 目 视 识别 撞击 坑 , 进而 到 近年 来 通过 计算 
1 自动 识别 撞击 坑 。 随 着 月 球 影 像 与 地 形 数据 的 大 量 获取 , 促使 撞击 坑 的 自动 识别 提取 成 为 
一 个 新 的 研究 热点 ， 近 年 来 ,众多 国内 外 学 者 提出 了 大 量 研究 方法 ， 并 建立 了 多 种 尺度 的 撞 
击 坑 数 据 库 中 [344516h1。 通 过 对 现 有 方法 的 分 析 总 结 ， 撞 击 坑 识 别 的 方法 大 致 可 以 归纳 为 
4 类: 人 工 目 视 识 别 、 基 于 形态 拟 合 、 基 于 机 器 学 习 以 及 基于 地 学 信息 分 析 的 方法 。 人 工 目 
视 识 别 法 : 利用 影像 和 地 形 数据 ， 建 立 识别 撞击 坑 的 解 译 标识 ， 在 计算 机 环境 下 ， 采 用 手工 
交互 的 方式 对 每 个 坑 逐 个 识别 , 这 种 方式 是 目前 最 准确 的 撞击 坑 识别 方法 。 这 种 方法 的 优点 
是 准确 性 较 高 ,但 是 费时 费力 ， 不 适合 于 全 月 球 和 小 尺度 撞击 坑 识别 。 基于 形态 拟 合法 : 由 
于 撞击 坑 在 图 像 上 有 类 圆 形 的 几何 特点 , 在 图 像 分 割 的 基础 上 以 霍 夫 变 换 和 二 次 曲线 为 代表 
的 形状 拟 合算 法 用 于 模拟 提取 完成 的 撞击 坑 边 界 ， 如 文 [8] 利 用 组 合 霍 夫 变 换算 法 ， 文 [9] 采 
用 模糊 霍 夫 变 换算 法 对 撞击 坑 进 行 识别 ， 文 [10] 提 出 一 种 基于 弱 中 点 霍 夫 变换 算法 ， 文 [1]] 
采用 最 小 二 乘法 拟 合 边 缘 椭 圆 的 方法 实现 对 撞击 坑 的 提取 。 基 于 机 器 学 习 的 方法 : 文 [12] 采 
用 BP 神经 网 络 加 训练 样本 的 算法 ， 文 [13] 利 用 决策 树 算法 ， 文 [14] 利 用 面向 对 象 方法 实现 
对 撞击 坑 的 识别 。 | 
等 利用 500m 分 辨 率 的 嫦娥 一 号 数字 高 程 模型 数据 和 形态 拟 合 的 方法 , 进行 撞击 坑 的 识别 提 
5 KUSARI, SOME ATTA, REEERE 

了 分 析 。 

总 结 前 人 研究 方法 的 优 缺 点 和 成 果 的 不 足 , 本 文 使 用 目前 世界 上 履 盖 最 全 、 精 度 最 高 的 
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嫦娥 二 号 数字 地 形 数据 , 提出 基于 等 高 线 地 形 分 析 的 方法 , 针对 小 型 撞击 坑 进 行 识别 提取 研 
究 ， 填 补 现 有 撞击 坑 数 据 库 小 型 撞击 坑 数 量 不 足 和 准确 率 不 高 的 缺陷 。 


1. 数据 


2008 年 2 H, 


国务 院 批准 探 月 工程 二 期 立项 ， 二 期 工程 主要 目标 是 实现 月 面 软 着 陆 ， 


开展 月 面 就 位 探测 与 自动 巡视 勘察 。 嫦 


二 号 作为 二 期 工程 的 技术 先导 星 , 首要 的 科学 任务 


是 对 全 月 球 和 嫦娥 三 号 的 月 面 着 陆 区 进行 详 查 , 精细 测绘 全 月 球 的 地 形 地 貌 。 嫦娥 二 号 卫星 
于 北京 时 间 2010 年 10 月 1 日 18 时 59 分 57 秒 在 西昌 卫星 发 射 中 心 成 功 发 射 。10 月 24 日 


16 时 49 分 ， 嫦 娥 二 号 了 


后 视 丙 个 视角 的 CBE o 


D 星 上 的 两 线 阵 CCD 立体 相机 首次 开机 工作 ， 并 成 功 获取 了 前 视 与 
至 到 2011 年 5 月 20 日， 嫦娥 二 号 卫星 CCD 立体 相机 完成 全 


部 任务 ， 共 获取 607 轨 图 像 数 据 。 全 月 球 影像 图 制作 过 程 中 选取 的 100km 轨道 高 度 成 像 的 
南北 纬 70° 范 围 344 轨 ， 极 区 共 40 轨 ) ， 用 于 全 月 球 数字 高 程 模型 数据 
E 流 程 主要 内 容 包 括 : 数据 准备 、 测 区 数据 组 织 、 平 差点 准备 、 全 球 平 差 、 


384 轨 数 据 (其 中 ， 
产品 生产 。 数 据 处 至 


测 区 初级 产品 生产 、 


据 、5 个 月 面 绝 对 控制 点 等 


50m 分 辨 率 数字 


标准 产品 生产 、 全 球 数据 灸 


高 程 模型 数据 。 嫦娥 二 号 数字 


整 性 镶嵌 精 度 等 方面 优 于 


据 进行 小 型 撞击 坑 的 识别 提取 。 


2. 方法 


2.1. 撞击 坑 模型 分 析 


根据 基本 形态 


E 棋 等 。 结 合 剥 娥 二 号 卫星 轨道 与 姿态 测量 数 


， 进 行 全 球 平 差 .经 过 4 年 的 数据 处 理 ， 成 功 制作 了 7m、20m、 


高 程 模 型 数据 在 空间 分 辨 率 、 数据 一 致 性 和 完 


国际 同类 全 月 球 数 字 产 品 。 本 文选 择 7m 分 辩 率 的 数字 高 程 模型 数 


， 撞 击 坑 可 分 为 简单 撞击 坑 〈 碗 型 小 型 撞击 坑 )、 复 杂 撞 击 坑 〈 具 中 央 峰 


的 较 大 型 撞击 坑 )、 多 环 撞击 坑 〈 撞 击 倪 地 )〉3 类 9。 一 个 理想 状态 下 的 复杂 撞击 坑 模 型 主 


要 由 坑 唇 边缘 (rh)、 坑 壁 (w)、 坑 底 ( 人 站、 撞击 培 
撞击 溅 射 堆积 物 由 坑 层 边 缘 (rh)、 丘 状 沉积 物 (rn) 和 辐射 线 沉 积 物 (rec) 构成。 撞击 溶 岩 (m) 获 盖 
在 其 他 形态 要 素 之 上 。 简 单 碗 型 撞击 坑 〈c) 一 般 分 布 在 丘 状 沉积 物 (rr) 之 上 。 撞 击 坑 识 


岩 (m) 和 中 央 峰 或 中 央 环 (p) 5 种 形态 要 素 组 成 。 


别 是 一 种 特征 识别 ， 对 影像 或 地 形 数 据 进 行 分 析 处 理 ， 识 别 撞击 坑 的 坑 层 边缘 (th) ， 并 提 


取 撞 击 坑 的 直径 CD) 、 中 心 位 置 和 高 程 等 信息 。 
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Fig. 1 Ideal complex impact crater model 

等 高 线 是 地 形 图 上 高 程 相等 的 相 邻 各 点 所 连 成 的 闭合 曲线 是 最 为 常用 的 表达 地 形 形 态 
特征 和 起 伏 情况 的 形式 。 等 高 线 可 以 从 数字 高 程 模型 中 提取 出 来 ,能 够 科学 地 反映 月 表 的 高 
程 、 坡 度 、 坡 向 、 高 地 走势 等 基本 月 表 地 形 地 瑶 。 嫦 娥 二 号 数字 高 程 模型 数据 的 空间 分 辩 率 


6 坑 模 型 (Don wilhelms, 1987) 了 2 简单 撞击 坑 和 复杂 撞击 坑 剖 面 图 


Fig.2 Impact crater structure 
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很 高 ， 月 表 高 程 在 二 维 地 理 空间 上 得 到 了 连续 表达 ， 利 用 该 数字 高 程 模型 数据 可 较为 方便 、 
精细 地 提取 月 球 撞击 坑 地 形 等 高 线 ， 进 而 表达 月 球 撞击 坑 形 态 ， 并 对 其 进行 分 析 和 处 理 。 


图 3 等 高 线 示意 图 
Fig.3 Contour diagram 

从 等 高 线 的 形态 分 析 ， 撞 击 坑 具 有 以 下 典型 空间 特点 (1) 撞击 坑 坑 层 边 缘 等 高 线 通 
常 为 近似 圆 形 或 椭圆 形 ; (2) 撞击 坑 层 坑 边缘 等 高 线 内 侧 通 常 存在 多 条 肉 套 的 相似 等 高 线 ， 
多 为 同心 圆 或 椭圆 ， 且 间距 较 小 ， 外 侧 等 高 线形 状 不 定 ， 且 间距 较 远 ,与 坑 层 边缘 有 明显 区 
别 ;，(3) 撞 击 坑 坑 底 高 程 小 ， 坑 边缘 高 程 大 ， 其 剖面 线 可 以 用 下 凹 的 抛物 线形 态 拟 合 。 因 此 ， 
通过 撞击 坑 地 形 等 高 线 的 特征 以 及 与 其 他 等 高 线 之 间 的 相互 关系 , 可 以 分 析 识 别 出 月 表 的 撞 
击 坑 。 从 图 4 的 撞击 坑 剖 面 形态 看 ， 等 高 线 表 达 的 是 撞击 坑 在 空间 上 的 等 值 形态 。 因 此 ， 只 
要 选取 合适 的 等 高 距 , 撞击 坑 唇 坑 边 缘 的 闭合 近 圆 形 等 高 线 就 可 以 盘 近 真实 边缘 ， 可 以 准确 
识别 提取 撞击 坑 的 地 理 位 置 、 大 小 和 深度 等 信息 。 


图 4 撞击 坑 等 高 线 示 意图 


Fig. 4 Crater contour diagram 


2.2. 指标 体系 和 参数 的 确定 


撞击 坑 等 高 线 往往 表现 为 相互 拒 套 的 类 圆 环 状 , 因此 , 初步 判别 撞击 坑 特 征 的 思路 是 根 
等 高 线 的 圆 度 、 和 矩形 度 、 体 态 比 等 指标 参数 保留 满足 条 件 的 等 高 线 , 再 分 析 等 高 线 之 间 的 
关系 ， 聚 合欢 套 的 等 高 线 并 标示 外 层 等 高 线 ; 然后 ， 通 过 剖面 线 的 凹凸 剔除 穹 隆 等 正 向 地 形 
区 域 ， 保 留 负 向 地 形 的 撞击 坑 区 域 ; 最 后 通过 霍 夫 变换 拟 合 撞击 坑 的 边缘 ， 得 到 撞击 坑 的 中 
心 位 置 和 大 小 。 因 此 ， 指 标 体系 和 参数 包括 圆 度 、 拖 形 度 、 体 态 比 、 谢 面 线 凹 凸 度 。 

() 圆 形 度 是 图 像 处 理 中 非常 重要 的 概念 ， 用 于 特征 提取 与 描述 ， 其 值 反映 了 被 测量 边 
界 的 复杂 程度 ， 公 式 为 


e=47A,/P* (1) 
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其 中 ，Ao 为 面积 ; P 为 周 长 ，e 为 44 面 积 与 周 长 的 平方 的 比 。 这 个 特征 对 圆 形 面 域 取 最 大 值 1， 
坑 层 边 缘 越 接近 1 的 等 高 线 越 符合 条 件 。 
(2) 矩形 度 是 常用 的 目标 物体 矩形 拟 合 因 子 ， 其 公式 为 


e=A,lA, > (2) 


JEP, AoA: ARAB IRIE ATTA: e ZEA FRAT RB) Sb TE HY FE ETE o 
对 于 矩形 物体 ，e7 取 得 最 大 值 1， 对 于 圆 形 物体 ，e 取 值 为 0.79; 对 于 其 他 细 长 物体 ，e 取 值 
变 小 。 甜 形 拟 合 因子 的 值 限定 在 0 与 1 之 间 。 

(3) 体态 比 定义 为 


e,=W/L ’ (3) 


RF, WA) eM RATE BE LAB) IMR FIBA » FALE UI AS LEH AL, RTE 
BUT DARE AS AS LU Aes CA BR o 

(4) HMRC: BAA eek ee eS, I wet 产生 近 圆 形 的 闭合 
等 高 线 。 潜在 撞击 坑 中 还 存在 , 对 此 情况 的 处 理 方法 是 判断 潜在 撞击 坑 所 在 的 数字 高 程 模型 
HRERS PUES, WRN FUES, WERTH; 反之 ， 若 剖面 线 为 -相形 态 ， 则 不 
是 撞击 坑 。 根 据 内 等 高 线 与 等 高 线 的 数值 大 小 关系 确定 是 否 为 真实 撞击 坑 。 

只 有 满足 指标 的 等 高 线 才 会 保留 下 来 ， 并 取 满 足 大 于 等 于 3 条 包含 关系 的 等 高 线 闭 合 等 
高 成 组 的 最 外 人 等 高 线 作为 撞击 坑 的 最 终 边 界 。 


2.3. 撞击 坑 边 界 的 拟 合 


对 满足 条 件 ， 被 保留 下 来 的 等 高 线 ， 利用 霍 夫 变 化 的 方法 拟 合 其 边界 ,最终 得 到 中 心 点 
和 半径 大 小 。 
霍 夫 变换 算法 是 图 像 处 理 中 识别 几何 形 以 及 其 特征 的 一 种 常用 方法 。 霍 夫 变 换 的 思想 
是 : 对 二 值 图像 T 上 任意 一 个 边缘 点 1(x,y)， 如 果 它 在 以 (a,b) 为 圆心 , 7 为 半径 的 圆 上 ， 则 满足 : 


2 2 _ 2 
(x-a) +(y-b) =r ， (4) 


从 (4) 式 可 以 看 出 ， 每 个 边缘 点 I(x,y) 映 射 成 空间 (a,b,r) 上 一 个 二 次 曲面 。 将 (4) 式 写成 
参数 形式 : 


a=x—rcosé | b=y-rsinO  ， (5) 


sep, 910,27) 。 设 待 检测 圆 的 半径 7s [名 ,名 ]， 则 按照 霍 夫 变换 ， 对 任意 一 个 边缘 点 


Kx y) 0 raa ,Ar 的 步 长 ， 遍 历 [0,2Z a IRI R21 ， 得 到 参数 空间 (a，b，n 的 


子 空间 ， 并 对 相应 的 累加 器 单元 AC (a，b，7) 加 1。 对 图 像 上 所 有 的 点 遍历 之 后 ， 得 到 累加 
器 数组 AC， 其 中 任意 一 点 4C (a, bs 表示 在 圆 (a,，b，7) 上 的 边缘 点 的 数目 。 相 同 数 目的 边 
缘 点 在 不 同 半径 的 圆 上 被 认为 是 圆 的 可 能 性 不 同 ， 半 径 越 小 可 能 性 越 大 。 因 此 定义 圆 度 的 
概念 : 


p=n/r ， (6) 


其 中 ，7 为 圆 的 边缘 点 的 数目 ，7 为 圆 的 半径 ，z 为 圆 度 。 将 累加 器 数组 4C 中 的 元 素 按 式 计算 
圆 度 ， 得 到 三 维 圆 度数 组 C。 圆 形 目标 对 应 于 圆 度数 组 的 局 部 峰值 ， 即 对 应 的 Co， ， 刀 为 
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3.1. 识别 流程 

在 嫦娥 二 号 7m 全 月 数字 高 程 模型 数据 中 随机 选择 了 G227, G228 区 进行 实验 , 基于 ENVI 
和 IDL， 按 照章 节 2 中 的 原理 和 方法 ， 编 写 程序 进行 识别 提取 。 首 先 ， 设 置 合理 的 等 高 距 在 
嫦娥 二 号 数字 高 程 模型 数据 中 提取 等 高 线 ; 然后 ， 通 过 圆 形 度 、 珑 形 度 、 体 态 比 、 剖 面目 凸 
性 ， si ei ie EOC, 最 后 ， 通 过 霍 夫 变换 提取 撞击 坑 边 界 信息 。 将 提取 的 信息 导入 
AirGIS， 确 定 撞击 坑 的 地 理 位 置 和 半径 大 小 ， 并 通过 高 程 差 计 算得 到 撞击 坑 的 深度 信息 ， 生 
成 具有 位 置 、 大 小 、 深 度 的 小 型 撞击 坑 数据 库 。 


FE Sone eae OO 


tras 


图 5 IDL 程 序 计 算 演示 


Fig. 5 IDL program diagram 
3.2. 识别 结果 
基于 数字 高 程 模型 数据 等 高 线 自动 识别 结果 如 图 6， 具 体 统计 结果 见 表 1。 


图 6 撞击 坑 识别 结果 示意 
Fig. 6 the crater detection result diagram 


本 文 的 方法 在 G227、G228 区 分 别 识别 撞击 坑 728 和 989 个 。 为 了 评价 本 文 方法 的 精度 ， 


使 用 Shufelt 8 提出 的 算法 评价 因子 : 检测 百分比 为 (Detection percentage)D=100TP/(TP+FN), 
分 支 系 数 为 B=FP/TP, 质 量 百分比 为 GO=100TPATP+FP+FN)。TP 表 示 正 确 识别 的 撞击 坑 数 量 ; 
oo 1 的 撞击 坑 数 量 ，FP 表 示 错 误 识别 的 撞击 坑 数 量 。 目 前 最 好 的 评价 方法 ， 依 
是 采用 人 工 识别 的 结果 作为 对 比 参照 。 
表 1 撞击 坑 识 别 结果 统计 表 


Table 1 Crater detection result statistics 


TP FN FP D B Q 


G227 697 53 31 92.9% 4.4 89.2% 


G228 918 81 71 91.9% Ld 81% 


从 表 1 统 计 的 识别 结果 可 知 ， 本 文 方法 在 两 个 实验 区 的 评价 因子 D 分 别 为 92.9%、91.9%， 
均值 为 92.4%，B 分 别 为 4.4%、7.7%， 均 值 为 6.1%，0Q 分 别 为 89.2%、81%， 均 值 为 85.1%。 
从 评价 因子 可 以 看 出 ,本 文 的 方法 对 于 小 型 撞击 坑 有 非常 好 的 检测 结果 ,正确 识别 率 很 高 ， 
且 错 分 率 很 低 ， 存 在 小 部 分 遗漏 识别 情况 。 


3.3. 讨论 


本 文 提 出 来 一 种 基于 高 精度 数字 高 程 模型 数据 的 派生 数据 等 高 线 通过 计算 机 自动 识别 
提取 小 型 撞击 坑 的 方法 。 相 对 于 传统 的 逐 点 分 析 坡 度 指 数 ,提取 图 像 边缘 转换 成 对 地 形 线 的 
分 析 ， 可 以 有 效 避 免 地 形 皱 裙 和 碎 石 群 的 干扰 ,识别 流程 更 为 简单 高 效 ， 识 别 速度 也 较 快 ; 
可 以 设置 不 同 的 等 高 距 , 分 批 次 分 区 间 识 别 不 同 尺度 的 撞击 坑 ， 其 识别 正确 率 高 ,错误 识别 
率 低 , 也 得 益 于 嫦娥 二 号 卫星 的 高 精度 数据 , 遗漏 识别 率 比较 低 , 说 明 该 方法 是 可 行 有 效 的 。 
不 过 需要 注意 的 是 ， 对 于 部 分 退化 严重 的 撞击 坑 ， 坑 缘 破 坏 严 重 , 坑 壁 平缓， 等 高 线 不 能 明 
显 表 达 和 准确 勾画 出 撞击 坑 的 形 貌 轮廓 ， 无 法 被 准确 识别 ， 导 致 漏 分 情况 的 发 生 。 但 通过 对 
比 同位 置 的 影像 数据 ， 通 过 人 人 工 识别 分 析 ， 是 可 以 分 辨 的 。 因此, 今后 进一步 的 研究 方向 和 
重点 在 于 如 何 识别 退化 严重 的 撞击 坑 , 这 也 是 所 有 撞击 坑 识 别 方法 中 的 注 弱 环节 和 有 竺 解决 
的 问题 。 在 实际 应 用 中 , 可 用 该 方法 的 识别 结果 补充 现 有 月 球 撞击 坑 数 据 库 中 的 小 尺度 撞击 
坑 ， 丰 富 现 有 的 月 球 撞击 坑 数 据 库 。 
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Automatic small crater recognition using digital elevation model 


from Chang’E-2 by contour line 


Yuxuan Liu!?3, Chunlai Li 13 
(LNational Astronomical Observatories, CAS, Beijing 100012, China) 
(2.University of Chinese Academy of Sciences, Beijing 100049, China) 
(3.Key Laboratory of Lunar and Deep Space Exploration, Chinese Academy of Sciences Beijing 100012, China.) 
Abstract: Crater detection is a significant research focus in planetary mission. Most of 
algorithm are focus on big crater cause of the data resolution is not precision enough. So this 
article propose a method about small craters automatic recognition from digital elevation model of 
chang’e-2 data, which is the highest precision lunar global data so far. Firstly, we analyze the 
feature of little crater model to get the parameter and establish index system. Secondly, we get the 
contour of test data and extract the candidate. Thirdly, we use the Hough-Transform get the craters 
information, including the coordinate and diameter. According to the result of test district, we can 
get the conclusion that the method is useful. The accuracy of the detection is over 85%, and we 
can set different parameter for different range of craters. The method and the result are good 
supplement to recent crater detection research, especially to the little size craters. It is the 
foundation of global carter database in the future work. 
Key words: Crater recognition; Chang’E-2 DEM; topography analysis; contour 
line 
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